K-inflationary Power Spectra in the Uniform Approximation 
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The advent of explicit Dirac-Born-Infeld (DBI) inflationary models within string theory has 
drawn renewed interest to the cosmological role of unusual scalar field dynamics, usually referred to 
as k-inflation. In this situation, the standard method used to determine the behavior of cosmological 
perturbations breaks down. We present a generic method, based on the uniform approximation, to 
analytically derive the power spectra of scalar and tensor perturbations. For this purpose, a simple 
hierarchy of parameters, related to the sound speed of the cosmological fluctuations and its successive 
derivatives, is introduced in a k-inflation analogue of the Hubble flow functions. The scalar spectral 
index and its running are obtained up to next to next to leading order for all k-inflationary models. 
This result relies on the existence of a well-motivated initial state, which is not trivial in the present 
context: having the wavelength of the Fourier mode smaller than the sonic horizon is indeed not 
enough and some conditions on the dynamics of the sound speed are also required. Our method 
is then applied to various models encountered in the literature. After deriving a generic slow-roll 
trajectory valid for any DBI model, simple formulae for the cosmological observables are obtained. 
In particular, the running, as the spectral index, for the so-called UV and IR brane inflationary 
models is found to be uniquely determined by the 't Hooft coupling. Finally, the accuracy of 
these cosmological predictions is assessed by comparing the analytical approximations with exact 
numerical integrations. 

PACS numbers: 98.80.Cq, 98.70.Vc 
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I. INTRODUCTION 

The primordial matter perturbations held responsible 
for growth and formation of large scale structure are com- 
monly traced back to quantum fluctuations of a scalar 
field if which should have dominated the energy density 
in the Universe at early times. The primordial power 
spectrum of both the scalar and tensor perturbations is 
a calculational output of the inflationary scenario, and 
the recent Cosmic Microwave Background (CMB) exper- 
iments have gathered considerable evidence in favor of 

it aaaiia. 

A compelling virtue of the inflationary paradigm is 
the fact that it can be sustained by a whole class of 
scalar field potentials V(ip), provided these exhibit char- 
acteristics (i.e. in their slope and curvature) in agree- 
ment with the "slow-roll" conditions: while inflation is 
under way, the potential must dominate over the ki- 
netic energy <f 2 /2, that is, InV should be flat enough 
not to accelerate the field quickly. The slow-roll regime 
for inflation is described by a hierarchy of parameters 
€i, assumed to be small, in which the primordial power 
spectra can be analytically expressed as a Taylor expan- 
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A large body of literature is devoted to the search for 
inflationary scenarios, in the sense that they should be 
naturally motivated by a high energy physics theory [13J ] . 
For example, the toolkit of ten-dimensional super-string 
theory and its various compactifications to an effective 
four-dimensional field theory has been used to design can- 
didate potentials [l4l Hoi. [r| . 

An interesting approach studies the effect of modifica- 
tions to the kinetic term of the inflaton (k-inflation) [17| ; 
in the perturbation treatment, these modifications mani- 
fest themselves as a (possibly time-dependent) "speed of 
sound" c s 7^ 1 (the speed of light being c = 1) for the 
scalar Fourier modes [la |. It turns out that string in- 
flation models where the inflaton field is an open string 
mode are typically of this kind, with the speed of sound 
being a function of the background geometry [ijj, [2CJ |. 
Hence inflationary model building in string theory com- 
bines both candidate potentials and non-canonical field 
evolution. 

The prime example of these scenarios are the so- 
called brane inflation models, where the inflaton f cor- 
responds to the position of a _D-brane within in a higher- 
dimensional manifold [2l|, [22|, [23| • Being an open string 
mode, if has a Dirac-Born-Infeld (DBI) action, which is 
essentially the square root of the induced metric on the 
brane. This metric, in turn, contains information about 
the chosen background compactification of the extra- 
dimensions through the position-dependent brane ten- 



sion T(<p). The inflaton potential V(ip) can be of various 
shape, and its exact calculation remains a matter of ac- 
tive research [24], [25j ■ Typically, it receives contributions 
involving T{ip), but may also be affected by finer geo- 
metric detail such as the presence of other branes in the 
extra-dimensional background geometry. As a natural 
consequence of the non-canonical interplay between po- 
tential and dynamics in the k-inflation case, one can no 
longer trust the intuition that flat potentials support in- 
flation. Through T((p), the warping of the background 
acts as a break on the field, allowing potential energy 
domination and accelerated expansion of the Universe 
even when In V is steep. 

As it is clear from the above description, the standard 
methods of slow-roll inflation break down in the case of k- 
inflation. In particular, the formula expressing the back- 
ground trajectory is modified since the shape of T(ip) 
(and not only the shape of the potential as in the ordinary 
case) affects the motion of the mobile brane. Another 
new feature of k-inflation, which is of prime concern for 
this article, is that one can no longer calculate the cosmo- 
logical perturbations' power spectra using the standard 
techniques. Indeed, as already mentioned above, the 
scalar perturbations have now a time-dependent speed of 
propagation which prevents us to integrate the equations 
of motion in terms of Bessel functions. Moreover, con- 
cerning the perturbations' evolution, the "usual" Hubble 
flow functions e^ do no longer provide a sufficient descrip- 
tion. Therefore, the main goal of this paper is to put 
forward a general formalism for k-inflation, resembling 
as far as possible the usual slow-roll formalism, where all 
these issues can be addressed in a consistent and unified 
way. 

The paper is organized as follows. In Sec. [Til we use a 
combined hierarchy (ej,<Jj) of Hubble and "sound" flow 
functions such that ej, Si -C 1, i > 1 is the analogue of the 
standard slow-roll approximation [20, [2y, [27], [28| ■ Then, 
using the background equations of motion, the slow-roll 
trajectory is expressed as a quadrature. Our new for- 
mula can be applied to any DBI model (but not to k- 
inflationary models in general) characterized by the func- 
tions T((p) and V(ip) and is valid under a single assump- 
tion, namely ei <C 1. In a next step, we carry on through 
the calculation of the k-inflationary perturbation spectra 
in full generality (but assuming, as usual, the smallncss 
of the Hubble and sound flow parameters), including the 
non-trivial effects induced by a varying sound speed, in 
the so-called uniform approximation [23, [3(| ■ In partic- 
ular, we derive the scalar spectral index and the running 
at the next to next to leading order for a general model 
of k-inflation. Moreover, we show that a time-dependent 
sound horizon may lead to sub-sonic Fourier modes start- 
ing their evolution out of the Wentzel-Kramers-Brillouin 
(WKB) regime. This generically results in oscillations in 
the primordial scalar power spectrum. At the end of this 
section, we also discuss how the new parameters can be 
used to generalize to k-inflation the classification of (sin- 
gle field) inflationary models described in Ref. [l2j ■ Then, 



in Sec. IIII1 our approach is applied to various exam- 
ple models encountered in the literature. In particular, 
we recover or extend previous results concerning power- 
law DBI inflation J31], [32|, the Kachru-Kallosh-Linde- 
Maldacena-McAllister-Trivedi (KKLMMT) model 111, 
chaotic Klebanov-Strassler (CKS) inflation [H, GH, H|] , 
and derive new results for models having a CKS potential 
plus constant |35|. At various stages, we assess the accu- 
racy of our approximation scheme from exact numerical 
integrations. In Sec. IIV[ we recap and discuss our main 
findings. Finally, in three short appendices, we briefly 
compare our new hierarchy of parameters to the param- 
eters already considered before in the literature and we 
recall how the master equation for the Mukhanov-Sasaki 
variable can be derived in k-inflation. 



II. DBI SLOW-ROLL INFLATION 

In this section, although we are concerned with the 
DBI models, the definition of the Hubble and sound flow 
parameters is valid in full generality for all k-inflationary 
models 17] . The same remark holds for the calculation of 
the primordial power spectra and will be made apparent 
by using explicitly the sound speed c s in the calculations. 



A. Basic equations 

The basic construction behind brane inflation models 
in super-string theory is discussed in several recent re- 
views [1J, [[(I [3g, |37| . Our starting point is the effective 
four-dimensional inflaton action 



S=- j ' d'x^—g 



1 



T (<p)\l 1 + 7 F 7^9^d fl <pd u i P 



V(<p) - T{fp) 



(1) 



with V(ip) the potential and T((p) the warp function. The 
warp function is determined by a specific choice for the 
geometry of the extra-dimensions. The shape of the po- 
tential receives many different contributions, which in- 
clude Coulomb-like terms describing the attraction be- 
tween branes and anti-branes as well as terms that arise 
from the embedding of different dimensional branes. On 
the string theory side, there is an ongoing debate on the 
number and form of these contributions. The form of 
T((p) is the subject of less controversy since known string 
inflation models use the singular conifold or its cousin, 
the deformed conifold [33|. However, T{ip) has also been 
treated as a completely general function in the literature, 
and here, both V(ip) and T(<p) will be considered as free 
functions in the sake of generality. 

It is worth noticing that the action (jTJ) defines a consis- 
tent theory. Indeed, for an arbitrary model of k-inflation, 



the action of the scalar field can be written as 

s k = [ d 4 xV=gP(X,v), (2) 

One can 



where the quantity X = — (\/2)g^d^pd v p. 
show 38] that this theory is well-defined if 



dP 



>0, 



2X 



d 2 P dP 



>0. 



(3) 



dX ' "' "" dX 2 dX 
The first condition comes from the requirement that the 
Hamiltonian should be bounded from below while the 
second is necessary if one wants the field equations to 
remain hyperbolic [38|]. In the case of DBI, one has P = 
—TaJI — 2X/T and one can check that both conditions 
are indeed satisfied regardless of the sign of the brane 
tension T(tp). In the following, we will always restrict 
ourselves to models of k-inflation that fulfill the above- 
mentioned conditions. 

Variation of Eq. (|TJ) with respect to the metric gives 
the DBI stress-energy tensor, 

2 5S 



T„ 



which is found to read 



T, 
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V(<p)+T(<p)[~-1 



where the Lorentz factor is defined as [2J 
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The relativistic analogy becomes evident in the case 
where one considers a spatially homogeneous field p(t) in 
a Friedmann-Lemaitre-Robertson- Walker (FLRW) uni- 
verse so that Eq. © simplifies to 



1 



(7) 



v/1 - p 2 /T( V ) ' 

where a dot denotes a derivative with respect to cosmic 
time. Clearly, y/T(p) here plays the role of an upper 
limit on the inflaton's velocity p. When expanding 7 for 
ip 2 <C T(ip), the action at first order resumes its canonical 
form. Moreover, the energy density and pressure read 



p=(7-l)Tfo>) + Vfa), p = 



7 



1 



7 



T(<p)-V(<p), (8) 



so that we have the Friedmann-Lemaitre equations 

K 



H 1 



[(7 



2H + 3H A = k 



1)T + V], 
1-7 



7 



T + V 



(9) 

(10) 



while the Klein-Gordon equation for the field reads 



.. 3H . 

p+^-p- 
T 



37 -7 3 



2 7 3 



2 AT 
dp 



1 dV 



7 3 dtp 



0. 



The constant k is defined by k = Sir /m 2 
four-dimensional Planck mass. 



(11) 
m pl being the 



B. DBI slow-roll trajectory 

In standard inflation, one usually defines a hierarchy 
of Hubble flow parameters from [I(J, [HI, fl2j 



Cri+l 



dln|e„| 
diV 



ro 



Hi 

H 



^T- (12) 



Their physical interpretation is that the expansion is ac- 
celerated as long as ei < 1 (potential energy domination). 
The slow-roll approximation assumes moreover that one 
has |e,| <C 1, i > 1, a condition which is in general nec- 
essary in order to have a sufficient number of e- folds. In 
addition, this last condition also allows us to integrate 
analytically the field trajectory and to compute the cos- 
mological perturbations' power spectra. 

In DBI inflation, we still retain the definition (fT5|) , the 
only subtlety being that the Hubble parameter is now 
given by Eqs. ([9]) and (fTO]) . Expressed in terms of deriva- 
tives of H with respect to ip, the first two Hubble flow 
functions read 



2 1 

kj H 2 



(2 



2 
«7 



dH 

dp 

2 (dH 
H 2 \~dlp~ 



2 d 2 H 
HJ^ 2 



1 d7 1 dH 

7 d(/? H dp 



(13) 
(14) 



In comparison with the standard case, we see that the 
expression of ei contains a 7 factor in the denominator. 
This merely expresses the fact that, even if the potential 
is not flat, inflation may occur provided 7» 1. Let us 
also recall that the above definition is fact valid for any 
k-inflation model with the replacement 7 = l/c s . 

In fact, one does not need more to derive the slow-roll 
trajectory for any DBI model. Let us first notice that 
Eq. (O can be recast into 



dN 
dp 



7 2 ff 2 



(7 2 - 1)T 



(15) 



The Lorentz factor 7 can now be expressed exclusively 
in terms of H and p. Indeed, Eq. ([TO]) together with 
H = ip dH/dtp yields 



V 



2 dH 

K7 dp 



(16) 



which can be used to replace p in Eq. ([7]). Solving for 7 
leads to 



7 ^ = v i+ ^0 



(17) 



Therefore, despite the fact that 7 contains a <p factor, 
it can be viewed as a function of the inflaton field only. 
Moreover, as shown in Ref. [4Q|, Eq. © can be recast as 



3H 2 = 



kV 



27 



(18) 



3(7+1) 



(1 



Let us notice that, up to this point, all equations are 
exact. To proceed further, we use the slow-roll approxi- 
mation to simplify Eq. (|18p by assuming £1 < 1. In this 
limit, since 7 > 1, one has 



H 2 



1 



-kV. 



(19) 



The DBI analogue to the standard slow-roll trajectory 
is readily obtained by replacing 7 in Eq. (|15p from its 
expression (ITT)) and using Eq. (|T9"]) for the Hubble pa- 
rameter, 



N{<p) 




(20) 



As a result, only the knowledge of V and T is required 
to calculate the DBI slow-roll trajectory, just as knowing 
the potential accomplishes the same goal in the standard 
case. Let us mention again that the only assumption 
that goes into obtaining Eq. (|20|) is e\ <C 1. As we will 
see in the next section, additional approximations are 
nevertheless required at the perturbative level. 

The new degrees of freedom introduced by the warp 
function T(ip) sug gest to define an additional hierarchy 
of parameters [2q. l4l| . In fact, in the same way that the 
€i encode the Hubble parameter evolution, it is conve- 
nient to consider their equivalent in terms of the "sound 
horizon". Therefore, we define the 5i, the sound flow 
functions, in a way similar to the Hubble flow parame- 
ters, but starting with the sound speed c s : 



>n+l 



din \5„ 
AN 



So 



^ in 



(21) 



In the case of DBI inflation, one gets for the two first 
parameters 

2 1 d 7 1 dH 



Si = - 



K7 7 dip H dip 



(22) 
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«7 



2 d7 1 dH 

7 dip H dip 



d 2 j/dip 2 1 dH 
d"f/dip H dip 



1 fdHY 



1 d 2 H 

H~dup 2 



(23) 



The full hierarchy is therefore given by the combined set 
(ei,di). Let us notice that, as in standard inflation, it 
is also possible to introduce various sets of flow parame- 
ters and one could also have a set of potential and warp 
function based parameters. Such an alternative hierarchy 
(e v , e T ) is summarized in the Appendix fA] for complete- 
ness. Finally, the e^ and Si defined here correspond to a 
subset of the "brane inflation flow functions" previously 
defined in Refs. M, M,M. 



1. Equation of motion for the scalar modes 

The main gauge-invariant equations for the cosmologi- 
cal perturbations are reviewed in Appendix [Cj There, 
it was shown that the (Fourier amplitude of the) 
Mukhanov-Sasaki variable obeys the following equa- 
tion [El 



<+ c 2 k 



,2,„2 



1> k = 0, 



(24) 



where a prime denotes a derivative with respect to con- 
formal time and where the function z is given by the 
following expression: 



cup 
~H 



1\ 3/2 12 a 
— =\ V £ i- 



(25) 



As in the standard case, one obtains the equation of a 
parametric oscillator. However, there is an important 
twist: the sound speed c s is no longer equal to unity but 
is now a time-dependent quantity which prevents the use 
of standard techniques to find solutions of this equation. 
In the DBI case, it is given by c s = I/7 < 1 and can take 
very low values. 

As in the standard case, one can entirely express the 
effective potential in Eq. ([24]) in terms of the Hubble and 
sound flow functions, 



[a^l/c s ] 
(a,y/e[/c g 



= H 2 



2-ei 



:£2 



1 2 1 
I 62 2 £1£2 



1 



+ (3-ei + e 2 )£i +61 + 6x52 



:£2£3 



(26) 



As expected, there are additional terms proportional to 
the sound flow parameters. 

As usual, the final goal is to compute the two-point 
correlation function or, in Fourier space, the power spec- 
trum. For the scalar modes, its expression reads 



A- 3 



k 3 



n-^ICul = 4?r2 



9 

a z t\ 



(27) 



where (k = v^/ z is the comoving curvature perturbation. 
In order to estimate this quantity, one has to integrate 
the equation of motion (f!M| from a set of initial conditions 
over the background solution, given by Eq. (|20[) for the 
DBI models. We now turn to the question of the initial 
conditions. 



2. Initial conditions 



C. K-inflationary perturbations 



The effective time-dependent frequency of Eq. (I24|) is 
^iven by 



In this subsection, we now turn to the theory of cos- 
mological perturbations in k-inflation. 



to 2 (k,ri)=c 2 k' 2 
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FIG. 1: Scalar power spectrum as a function of the number 
of e-folds Nk at which a given mode k crossed the sound 
horizon in the case of the chaotic Klebanov-Strassler model 
(see Sect. IIH B 2|l . The oscillations on large scales come from 
the violation of the WKB condition while the initial state is 
enforced to be a Bunch-Davies vacuum. The bottom panel 
shows the sound flow parameters as a function of the number 
of e-folds: 5\ ~^> 1 initially triggers the WKB violation. 



A well-defined and well-motivated initial state can be 
chosen in the adiabatic regime for which a WKB solu- 
tion exists, i.e. for 



Q 



<i 



where the quantity Q is 



W 2 

Auj 2 



a 



Q( k ^) = i-2- — 



U! 

2w 



(29) 



(30) 



In the standard case, the modes of astrophysical inter- 
est today are, at the beginning of inflation, such that 
their wavelength is smaller than the Hubble radius. This 
implies that u> ~ k and the quantity Q vanishes. As a 
consequence, the condition \Q/lu 2 \ <^C 1 is obviously sat- 
isfied and the WKB state 



exp 



«k(?7) 



i / u!(k,T)dr 



JHv-mm) 



VMM) 



(31) 



is the preferred initial state. 

However, in the k-inflationary case, the time- 
dependence of the sound speed brings new complications 
into this question [18|, [2CJ, [23, |42j . Indeed, even if initially, 
the modes are within the sound horizon and such that 
u> ~ c s k, the effective frequency is still a time-dependent 
quantity and, therefore, it is a priori not obvious that 
a well-defined state can be chosen in this context. One 



has to check for each mode that the quantity Q/oj 2 is 
indeed small. From Eqs. (f2"5|) and (f5U|) . this one can be 
expressed in terms of the sound flow functions and reads 

" ;/ f h-e l 5 l +8 1 8 2 + \8i). (32) 



U) 



2c 2 . k 2 



Assuming the slow-roll conditions on the e^ are satisfied, 
the WKB condition (|29|) can still be violated as soon as 
c s k/TL ~ 8±, and thus even for the modes deep inside the 
sound horizon provided 5\ is big enough. As an example 
of such a situation, Fig. Q] shows the scalar power spec- 
trum obtained from a numerical integration of Eq. (|24|) 
in a case where 8\ 3> 1 initially, the warp function and 
potential being those of the chaotic Klebanov-Strassler 
model, see Sect. IIIIB 21 Enforcing the modes to start in 
the Bunch-Davies vacuum is no longer justified and leads 
to oscillations in the power spectrum. Clearly, such a sit- 
uation essentially concerns the modes crossing the sound 
horizon soon after the beginning of inflation, the ones 
for which c s k/7i cannot be chosen big enough to satisfy 
Eq. (221. 

The conclusion is that, in order to be able to choose a 
well-motivated initial state for a given Fourier mode in 
k-inflation, it is not enough to have a wavelength smaller 
than the sound horizon. Additional conditions on the 
sound flow parameters 8\ and 82 are also necessary. As 
shown before, a situation where A& = 2?:a/k is smaller 
than the sound horizon and Q/uj 2 3> 1 can easily been 
designed. In such a case, our ability to work with a well- 
defined initial state is lost. 



D. The uniform approximation 

At first order in the Hubble and sound flow functions, 
one can use the uniform approximation to solve the mode 
equation. The uniform approximation was developed in 
Refs. [23, |3CJ. For this purpose, it is convenient to re- 
write Eq. (gl]) as 

.2 



c 2 k 2 



1/4 



V 



vu = 0, 



(33) 



where the function v (jf) can be calculated from the effec- 
tive potential given in Eq. (|2"6")l . Following Refs. [2g.l3fj|. 
one also defines the two following functions: 

9(v)^^-clk\ q(v) = ~^- (34) 

The so-called turning point is defined by the condition 
g{v*) — an d, for each mode, occurs at the time ?7*(fe) 
such that 

fcry* = — . (35) 

c s* 

The uniform approximation tells us that the Mukhanov 
variable i>k can be expressed as [29|, [30] 



Vk(rj) = Ak I - 
\9 



1/4 



.Y«f)+B k [£ 
\9 



1/4 



Bi(/) , (36) 



where A^ and B^ are two constants to be determined 
from the initial conditions and Ai(x) and Bi(cc) denote 
the Airy functions of first and second kind respectively. 
The function f(k, 77) is defined by 



f(k,v) = 



\y-y*\ 
77-77* 



drv^Ml 



2/3 



(37) 



As already discussed above, we assume that adiabatic- 
ity is valid initially then choosing the initial conditions 
as an initial state of the WKB form yields 



A k = iB u 



B k = J-c 



(38) 



where 9 is just an unimportant scale-independent phase 
factor that will be ignored in the following. The solu- 
tion (j3"6")l is now completely specified. 

After the turning point, using the asymptotic behavior 
of the Airy functions, the solution can be written as 



■"'"■■ WTJ^O' 2 



a 



(39) 



and one still has to calculate the integral in Eq. (f3~7| . For 
77 > 77*, it simplifies to 



;/ 3/2 m = 



dr 



c?k 2 . 



(40) 



The functions ^(77) and 1/c s (t7) can be expanded in terms 
of the Hubble and sound flow functions. For instance, as- 
suming that l/c s admits a polynomial expansion around 
77*, one has 



1 



c s (t?) 



E 

n=0 



1 



(n) 



(77 - 77, ) r 



At first order, one finds from Eqs. (fl"2"|) and (|2ip 



(n) 



— (n-l)\H r ' 



0(eS) , 



(41) 



(42) 



where 0(e6) stands for all terms of order two in the e,, Si 
or mixed. Let us notice that all terms in Eq. (|4"Tj) should 
be considered at first order in e, and 5i . Plugging Eq. (j4"2"|) 
into Eq. (J4TJ) yields an infinite sum that can however be 
resummed into 



c s (t7) 



1 - 5i*ln[l -7Y*(77 - 77,)] 



(43) 



a star indicating that the corresponding quantity is evalu- 
ated at the turning point defined above. For consistency, 
W* has still to be expanded. Using 



H(v) 



1+ei 



G(eS) 



one finally gets 
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c s(v) 



0{eS) 



(44) 



(45) 



Along the same line of reasoning, one would show that 
the function 7/(77) reads 

v 2 {v) = |+3ei*+^e 2 *+3<5i*+0(e<r) = ^+0(eS) . (46) 

One can check that this expression matches with the 
standard result by setting Si* = (see Ref. jTij). Let 
us notice that one may avoid the infinite re-summation 
by performing an expansion directly in terms of the num- 
ber of e-folds. Indeed, a Taylor expansion of l/c a (N) is 
also a flow expansion. For instance, one has 



1 = — + —S u (N - JV») + 0{5 2 ) 



c s (N) 



(47) 



To recover the conformal time dependency one may use 
Eq. (gU) to get 

?7* 



N-N^ln^^ln. 1 -^ 1 '" 



, In ^+ 0(e). (48) 

Plugging the previous equation into Eq. (|47[) immediately 
gives Eq. (|4"5|) . 

Let us now return to the calculation of the func- 
tion f(k,rj). Inserting Eqs. (|45|) and (|46|) into the in- 
tegral (I40p . and defining the new variable 



c s *krj 



(49) 



one obtains at first order in the Hubble and sound flow 
parameters 



^~ 



In 



1 + VT 



w' 



-5i*v 1 — 7« 2 ln(— w). 



(50) 



where we have used Eq. (I35[) . This expression will be 
used in the following to derive the power spectra. 



E. Power spectra and spectral index 

1 . Scalar power spectrum 

We are now in a position to estimate the power spec- 
trum of scalar perturbations. The comoving curvature 
perturbation being constant on super-sonic length scales, 
see Appendix [C] the power spectrum is obtained from 
Eq. pf|). using Eqs. (J39J) and JSOJ) in the limit w -> 0. 
Lengthy but straightforward calculations give, at first or- 
der in the Hubble and sound flow functions, 
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H: 



7rm^ei*c s 
1 



(I8e" 3 ) 
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(51) 



Notice that the above expression is still an implicit func- 
tion of k through its dependence on 77,. Let us also re- 
mark the presence of the factor 18e -3 ~ 0.896 in the over- 
all amplitude, which is typical for the WKB and uniform 
approximations. As discussed in Ref. [43|, this is rather 
unfortunate since this factor damages the approximation 
of the overall amplitude down to the 10% level. However, 
as shown for instance in Ref. 44 1. this problem can be 
rather easily fixed. Very roughly, one can renormalize 
18e -3 to one to recover the exact amplitude [43j| . The 
spectral index is, on the contrary, predicted accurately 
by the WKB and uniform approximations. 

To remove the implicit dependence in k hidden in 77*, 
we define a pivot wavenumber fc and expand all terms 
around an unique conformal time 770 which is the time 
when the pivot scale crossed the "sound horizon" , namely 



koi]^ 



(52) 



Then, one can use the flow expansions for TL, l/c s , and 
the 6j, Si, but now around ?/ . In this case, Eq. (|51j) 
becomes 



T>t 
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7rm^eioC s o 

+ {D + 2)(5 i0 - (2ei<> + e 2 o - 610) In 



1 - 2(D + l)eio - De 2 

k 



A'-<> 



(53) 



where we have defined D = 1/3 — In 3. From the above 
expression, one can already read off the spectral index of 
scalar perturbations, 



1 — — 2ei — €20 



(54) 



The standard expression is corrected by a term, <5i<>, 
which takes into account the time-dependence of the 
sound speed. 



2. Tensor power spectrum 



which is obviously different than 7/0 in Eq. (|52p . It is 
therefore convenient to express H^ and ei> in terms of 
the parameters evaluated at 77 = ?y . Using Eq. (|48| with 
V> = c s«^; a t fast order in the sound flow parameters, 
one gets 



^•-! = ^?( 18c ~ 3 ) 



7T?Tl~ 



1 - 2[D + l + ln—)( U 



2ei<> In ■ 



(57) 



Let us remark that, in the amplitude of the power spec- 
trum, the coefficient in front of the first slow-roll parame- 
ter is different as compared to the standard case. Indeed, 
in the standard case, c s0 = 1 and ln(l/c s0 ) vanishes. 

Eqs. ([53")) and ((57)) constitute one of the main result of 
this article. They represent the general scalar and tensor 
power spectra of k-inflation, valid at first order in the 
Hubble and sound flow parameters. In Ref. J45j, they 
have recently been compared to the WMAP5 data [l|, 0, 

IIM]. 



F. Running of the spectral index and higher order 
corrections 

Having obtained the scalar perturbation power spec- 
trum from the uniform approximation, we will now put 
this result to use in calculating the spectral index n s and 
the running a s of the spectral index using the method 
of Ref. [12}. Expanding the scalar perturbation power 
spectrum around the pivot scale fc<> in terms of ln(/«/£; ), 
one has 



W) = n(ME?f ln " f 



n>0 



(58) 



At zeroth order, the spectrum is given by the leading 
term of Eq. (fS"3")) . namely 



Although the equations of motion and evolution of the 
tensor perturbations are not affected by a varying sound 
speed in the scalar sector, there is however a subtle ef- 
fect associated with the choice of the e-fold at which one 
evaluates the Hubble and sound flow functions entering 
the Taylor expansion of the power spectrum. 

As shown in Ref. [431] . around the pivot scale k , the 
tensor power spectrum in the WKB, or uniform approx- 
imation, at first order in slow-roll, reads 



V h (k) = ^$ (I8e~ 3 ) 



l-2(Z» + l)ei > -2ei > ln — 



(55) 
where all background quantities are evaluated at the time 
77c, such that 



k<>Tl[> 



1. 



(56) 
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(59) 



Since the physical power spectrum must not depend on 
the pivot scale, d7 ? < r(fc)/dlnfc = 0, one may establish 
the recursion relation |12l ] 

d In V^ do„ 

a n+l = —r-~ — : — a n 



d In k d In k 

1 ( da n dlnP c \ 

H TT7 — a„ , n > 0, 



1 - eio + <5 i0 \ dA^o dA 



(60) 



where in the last expression we have used 

d\nk = (l-e l0 + S l0 )dN . (61) 



Using Eq. (f59|) . one gets 



and it can be identified with 



din 7? 



dN 



= -2e 



lo 



£2o + $k 



(62) 



In terms of the expansion ([58"| . the first coefficient ao is 
determined by the spectral amplitude, while a\ is related 
to the spectral index, and a 2 to the running. Note that 
if we know ao to gth order in (e^,^), we can determine 
a n to the order q + n. From the uniform approximation 
[see Eq. ([53")) ] we know that 

a = 1 - 2{D + l)ei<> - De 20 + (D + 2)<S l0 , (63) 

at first order. The recurrence relation, up to second or- 
der, therefore gives 

<zi = -2eio-e 2 o + £io + 2(2£> + l)e? 

-(AD + 3) ei (5i + (2D - l)ei e 2o + De 2 2<> 

-De 2o e 3o - (2D + l)e 20 6 l0 + (D + 1)& 

+ 0D + 2)«5i o <5 2o . (64) 

One more iteration of Eq. (|60[) allows to determine a 2 up 
to third order, 

a-2 = Ae 2 1(> + 2ei e 2 o - 4ei oi + e 2o - 2e 2<> 5 1<> - e 2o £3o 
+S 2 1o + S 1o S 2o + 6e? o e 2o + 12De 2 l0 6 l0 
-(2D - l)eioe 2o + 3(2D - l)eio£ 2 oOio - 6Deio<5? 
+3De 2 2 Ji<> - 3De 2 o^ + Do^ - 3De 2 o£3o<Sio 
+3(£> + l)<0 2 o - SDef^ + 2(2D - l)e l0 e 2<> e 3(> 
-6(D + l)eio<5io<5 2 o - De 3 2o + 3De 2 2o e 3<> 



-3(D + l)e 2o <Sio<5 2o - De^e^ - De 2c ,e 3< ,e ic 
+ (D + 2)6io6^ + (D + 2)o 10 <5 2o <5 3o . 



(65) 



The relation between the a^ and n s is obtained using the 
definition 



lis - f = 



din 7? 
dlnfc 



(66) 






which, when compared to the corresponding derivative of 
Eq. ([58]), leads to 



1 = ^-. 

a 



(67) 



and at second order in the Hubble and sound flow pa- 
rameters, one gets 



n s - 1 = -2ei<> - e 2o + S l0 - 2f? Y<> - (2D + 3)ei<>e 2c 



-3ei (5io + e 2 o<5io — -De2«£3o — 6 l0 
-(D + 2)6x0620- 



(68) 



The same calculation can be repeated for the running. 
Its definition reads 



d 2 lnP c 
din k 



(69) 
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Up to second order, one finds 

oi s — — 2eio£2o ^ £20^30 + <5io<52o . 



(70) 



(71) 



We see that the parameter 62 appears in the above ex- 
pression. Using the tools developed previously, one can 
even estimate the running at third order. The result 
reads 



2ei e 2 o — £2o£3o + <^io<52o — 6e 10 e 2o + 5ei e 2o <5i 
(2D + 3)e l0 ei - 2(D + 2)e 10 e 2o e 3o 
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Finally, from the power spectrum of the tensor pertur- 
bations, the tensor to scalar ratio at first order in the 
Hubble and sound flow parameters reads 






— 16c q o£lo • 



(73) 






We recover that since c s < 1, r is reduced compared to 
the canonical single field dynamics. 

We can compare our approach to the existing results in 
the literature, and, at first order, our results agree with 
those of Refs. [271 |28l [34 . l35j |. However, we would like 
to stress that at higher order, and in particular for the 
running, our results match only with those of Ref. [271 ]. 
Indeed, in Refs. [H, [H, [H[ , Eq. ([24]) is solved in terms 
of Bessel functions along the lines of the standard for- 
malism which assumes that c s is a constant. This ends 
up being an acceptable assumption at zeroth order only 
as is clear from Eq. ([47)) . The spectral index at first or- 
der being, roughly speaking, the derivative of the power 
spectrum amplitude at zeroth order, one may indeed ob- 
tain its correct first order expression with the assumption 
that c s is constant. However, if one wants to derive its 
quadratic corrections, or the running, then it is neces- 
sary to know the correct amplitude of the power spectra 
at first order. In Ref. [231, this goal was achieved using a 
conveniently chosen transformation of the time variable 
to absorb the time dependence of the sound speed, the 
resulting equation being integrable at first order. 

This argument is at the heart of the present article. If 
the sound speed of the perturbations is a time-dependent 
quantity and if the first order expression (|47[) of c s is in- 
serted into Eq. ([M|) . then the solution cannot be found by 
the usual technique. Therefore, the standard approach 
cannot be used and this prompts the use of a different 
method. This is what was done in Ref. [271 ] . using a new 
time variable, and what is done in the present paper un- 
der the WKB/uniform approximation. 



G. Model classification 



III. EXAMPLE MODELS 



To conclude this section, we generalize the classifica- 
tion of inflationary models of Ref. [12| to the DBI case. 
The energy density of a DBI inflaton field of Eq. IJ5J) can 
be re-written as 



7 -2 

<p 

7+1 



VM 



(74) 



from which it is easy to see that we recover the standard 
expression in the limit 7 — > 1. Let us refer to the first 
term in Eq. (|T4[) as the kinetic energy contribution, while 
the second term represents the potential energy. The 
Hubble flow functions e± and 62 then can be used to study 
the respective evolution of these contributions. With 
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(75) 
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which represents the generalization of Eqs. (6) and (7) of 
Ref. [12| , the change in the potential energy is given by 



V 



-Hjip 2 



' (e 2 -2e x ) ' 
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;<h 



(77) 



Note that again the standard expressions are recovered 
for 7 — > 1 since, in this case, also 8\ — > 0. From Eq. ([77]) . 
we see that the potential energy density can never in- 
crease for small values of e, and Si even if 7 is large. 

The change for the ratio between kinetic and total en- 
ergy density is given by 



dt 



m 



= H h- 



(78) 



Therefore, 62 = is the borderline between the regime 
where the kinetic energy contribution to p increases 
(f2 > 0), and the regime where the kinetic energy con- 
tribution decreases (e2 < 0). For a refined classification 
of DBI inflationary models, let us also calculate the time 
derivative of the kinetic energy density in Eq. (f?4")) : 



7 .2 
<P 



dt V7 + 1 



^J—^H^[{ 1 +l){e 2 -2e l ) + 5 l ] 



This equation can be viewed as the equivalent of Eq. (10) 
of Ref. [12[ . Hence, the kinetic energy density increases 
while 



62 > 2ei 



Si 

7+1 



(80) 



and decreases otherwise. We see that the standard condi- 
tion, 62 — 2ei > or e-2 — 2e\ < is modified and that the 
sound flow parameter 8\ now participates in the new cri- 
terion. This is natural since the factor 7 appears in the 
expression of the kinetic energy. However, in the limit 
7 — > +00, the standard condition is also recovered. 



We now illustrate our results in the case of DBI infla- 
tion with some specific choices of V(<p) and T(ip) con- 
sidered in the literature. The DBI slow-roll trajectory 
permits to calculate the resulting values of 7 and the 
first ti and Si parameters, and therefore the shape of the 
scalar primordial power spectrum through Eq. (|53|) . As 
a warm up, we start with power-law inflation, where, as 
it is the case in the standard situation, exact solutions 
are available. Then we turn to the more important case 
of branc inflation that we discuss in some detail. 



A. DBI power-law inflation 

The DBI analogue of power-law inflation has been 
studied in Refs. [3l|, [32J. As we show in the follow- 
ing, it is particularly convenient for testing the previous 
approximations since the perturbation equations are ex- 
actly solvable in this case. 

Looking for a power law behavior of the scale factor 
in terms of the conformal time, one finds that the warp 
function and the potential are given by [3l|, [33, |46| 



T(tp) = T exp 
V(<p) = Vo exp 



2 7 k 

(P-Vo, 

P 



27K 
P 



(<P - <Po) 



where the two constants To and Vo are related by 

7-1 " 
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(81) 
(82) 

(83) 



The solution of the Einstein equations are such that the 
Lorentz factor 7 is a constant and the scale factor and 
scalar field can be expressed as 



a(t) 

In this case, one has also the relation 

2^7 




T n = 
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(84) 
(85) 
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In terms of conformal time 77, the scale factor is also of 
the power-law form 



a(v) = £0 \V 



1+/3 



(87) 



with 



2p- 1 

P-I 
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One can check that when p — > oo, goes to —2 and one 
recovers the de Sitter case. 

Let us now turn to the calculation of the perturbations. 
The equation of motion (j2~4"|) for the Mukhanov-Sasaki 
variable takes the form 



i 2 



0G?+i) 



?r 



vu = 0, 



(89) 



where 7 is now a constant. This equation can be solved 
explicitly in terms of Bessel functions and the solution 
reads 



Vk(v) = (kr]) 1/2 

+ -Bk</-(/3+l/2) 



AkJp+l/2 



(90) 



where A^ and i?k are two scale-dependent arbitrary con- 
stants. As usual, one requires the initial state, evaluated 
in the limit krj/'j — ► —00 to be of the WKB type. This 
completely fixes A^ and B^ which are given by 



A k = -B k c l ^ +1 ^ 



Bu = 



r(/3+l/2)/2-Mr/4 

cos (tt/3) 



(91) 
(92) 



The function ([90]) is now completely specified. Using the 
asymptotic behavior of the Bessel functions in the limit 
h]/-f — > 0, one obtains the power spectrum from Eq. |2"7|) 
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where the function h{/3) stands for 

TV 
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(94) 



In the above expression, T is the Euler function of first 
kind. In the de Sitter case, h(— 2) = 1. Let us also notice 
that, in the case of power-law inflation, the first slow-roll 
parameter t\ is actually constant and given by 



(1 



2 + P 
1 + /3 



(95) 



As usual, for j3 = —2, one recovers t\ = 0. 

We can now compare the previous expression with the 
one obtained in the uniform approximation. The inte- 
gral (|40]) reads 



where, as already mentioned above, 7 is, in the present 
context, a constant. This integral can be performed ex- 
actly. Then, straightforward calculations lead to 



n 



iM 



l&i V7 
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where the function j(j3) can be expressed as 

2e 2/3+l 



m 



(2/3 + l) 2 ^ 2 ' 
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The above power spectrum is the analogue of Eq. (58) in 
Ref. [43|. In particular, the function j(/3) is exactly the 
same as in the standard case. Therefore the error in the 
amplitude of the power spectrum in the WKB/uniform 
approximation is similar (see Fig. 1 of Ref. [43] ) . In par- 
ticular, and as already mentioned, it is around 10% for 
p ~ —2, i.e. close to scale invariance. For (3 < —2, when 
the spectrum is not necessarily close to the Harrisson- 
Zeldovitch spectrum, the accuracy of the WKB/uniform 
approximation becomes better. Let us also mention that 
Ref. [44J has discussed how to improve the precision on 
the amplitude and the method developed in that article 
may be applied here. As can be noticed on the above 
equations, it turns out that the spectral index is pre- 
dicted exactly in the case of (DBI) power-law inflation. 



B. Brane Inflation 

We now apply our formalism to the case of the brane 
inflationary models discussed in Refs. [42l |47J. lij. |49|. |5C| . 
In these scenarios, the warp function and potential are 
given by the following expressions 



TM-*-, 
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(99) 



Here, the factor e stands for e — ±1 and the positive sign 
corresponds to the so-called "Ultra- Violet" (UV) mod- 
els while the minus sign refers to the "Infra-Red" (IR) 
scenarios [42|. The Coulomb potential is due to the at- 
traction between a _D3-brane sitting at the bottom of the 
throat and a mobile D3- brane, for a review see Ref. [4fJ| . 
The quadratic correction in Eq. (|99[) has a different sta- 
tus. It is a phenomenological description of the brane 
moduli potential and its shape is not established as neatly 
as for the Coulomb potential. The above described model 
is characterized by four parameters, the mass m, the scale 
jjl, the dimensionless 't Hooft constant A and Vq (whose 
dimension is m 4 ). In the following it will become clear 
that the evolution of the field is in fact controlled by two 
dimensionless parameters a and p. The parameter a is 
defined by 



127T771 2 
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96tH 
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(100) 



while the other dimensionless parameter /3 is given by 

V 



P 



(101) 



Physically, /3 measures the importance of the constant 
term relative to the mass term in the potential. 
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For those models, the vacuum expectation value (vev) 
of the inflaton field possesses a geometrical interpreta- 
tion, namely the distance between the two branes living 
in the throat. As a consequence, it can be written as 

p = \fTj,r, where T 3 = 1/ (27r) g s a' 2 is the tension 

of the brane, g s being the string coupling, a' the string 
scale and r the radial coordinate along the throat. Notice 
that we do not consider a possible angular motion of the 
brane [5ll |52I | . In order for the moving brane to be inside 
the throat, one should impose r < r where r is the 

uv uv 

coordinate at which the throat is connected to the bulk. 
This quantity can be written as [53| 



1. KKLMMT model 



47r<? s a' 



M 



(102) 



where A/" is a positive integer representing the total 
Ramond-Ramond (RR) charge and v represents the di- 
mensionless ratio of the five-dimensional volume forming 
the basis of the six-dimensional conifold geometry to the 
volume of the unit five-sphere. Requiring the volume of 
the throat to be smaller than the total volume of the 
extra-dimensions amounts to 



<P< V 



< 



uv v / 2^A7' 



(103) 



and inflation always occurs for sub-Planckian values of ip. 
On the other hand, the bottom of the throat is located 
at ro and, hence, one must have ip > (pa = \fT^r^. More- 
over, for the model to be valid, the (physical) distance 
between the brane must be larger than the string length 
and one can show that this amounts to 

ip > V2 st rg = <p e aVuv ■ (104) 

Notice also that the parameters of the potential ([55]) 
can be calculated in terms of the stringy parameters. The 
tension T(ip), in the model under consideration, can also 
be written as T(tp) = Ts(p/(p ) 4 which implies that |42j | 



A = 



2tt 2 v 



The constant term Vq is given by Vq ■■ 
can also be expressed as 

V = 2hHr )T 3 , 



(105) 
47t 2 u</?q/A/' which 



(106) 



where h(ip) = ip/ip is the warp factor. On the other 
hand, the constant /i can be expressed as fi 4 = p\jM . 

Finally, from the requirement that the volume of the 
throat is smaller than the total volume of the extra- 
dimensions, one deduces that [431 




1 h 2 (r ) 



V2A^ M 
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(107) 



We now study the different types of inflation that can 
occur in the scenario under consideration. Each type 
corresponds to a different choice of the parameters char- 
acterizing the model as discussed in Ref. [42j where the 
overall situation is summarized on two phase diagrams. 



The quadratic term in Eq. (|99[) can be neglected in 
the case where Vo/m 2 ^$> m 2 > p 2 or, in other words, 
for /? > 1. This, in turn, implies that a 3> 1 because 
of Eq. GH7L Tni s limit corresponds to the KKLMMT 
model [19|, |4(|. The trajectory given by Eq. (j2"0)) reads 
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(108) 

The field rolls down the potential (i.e. ip is decreasing) 
while inflation is under way and we always have <p > \x. 
Therefore, the previous expression can be approximated 
by 
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(109) 



For p 3> /i, the first two terms in this expression clearly 
dominate. In fact, this is precisely the trajectory found 
in Ref. [4(| in the limit where the DBI dynamics can be 
ignored. Therefore, the last two terms can be understood 
as the DBI corrective terms induced by the non-standard 
kinetics. Clearly, these corrections do not play an impor- 
tant role in the case of KKLMMT as long as ei <§; 1. As 
shown in Ref. [4(| , the DBI effects can only be significant 
for ei > 1 and we therefore do not proceed further with 
the KKLMMT model 1 . 



2. Chaotic Klebanov-Strassler models 

We now move on to the model where the quadratic 
correction dominates (hence e = +1) the potential. This 
corresponds to the condition j3 -C 1 and to the scenario 
considered in Refs. [2fJ,[39(, which we refer to as chaotic 
Klebanov-Strassler (CKS) inflation. The potential and 
warp function in this case read 



V(<p) 



1 
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n<p) 



A 



(110) 



The two remaining parameters are the mass m and the 
dimensionless constant A. The integration of Eq. (J2"D|) 



1 It was also shown in that reference that, depending on the values 
of the string parameters a' and g a (i.e. the string coupling), the 
end of inflation can occur either by violation of the slow-roll 
conditions or by instability. In the first case, the spectral index 
is given by n s ~ 0.97 while it is n B ~ 1 in the other situation. 
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can be performed explicitly and gives 
a 1 ' 2 



2 TV 
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where we have defined 

x 4 = a((^/m pl ) 4 . 
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(112) 



This is the implicit slow roll trajectory; in general, invert- 
ing this expression to find ip(N) is analytically impossi- 
ble. However, tp decreases during inflation, therefore the 
initial field value <£;„; is necessarily larger than all the field 
values ip attained during inflation. Assuming ip -C pi n i, 
an approximate inversion gives 



ip ~ ipi n i exp 



,1/2 



N 
4tt 



(113) 



This expression is nothing but the solution ip — ► 1/t 
found in Ref. [39( but expressed in terms of the number 
of e-folds. 

We can also derive the behavior of 7, t% and Si. Using 
Eq. ([T7]) together with Eq. (JTSJ) in the limit where ei < 1, 
we find the following expression for the Lorentz factor: 



— (^) J& + 1 
ip ) y mi. a 



(114) 



Note that 7 — > 00 as <p — > 0, i.e. the ultra-relativistic 
limit is attained at late times. Likewise, with the help of 
Eq. (HID, we find for e x 
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Note that ei approaches a constant as 93 — » 
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We recover that, in order to have slow-roll inflation for 
small values of tp, the parameters of the model, m and A, 
must be such that a <C 1. 

The above expression has also an important conse- 
quence with respect to the end of inflation in this model. 
Indeed, the condition t\ = 1 reduces to 
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Since we are in the limit a<l, this means that infla- 
tion cannot stop by violation of the slow-roll conditions. 
It must be stopped by another mechanism, typically by 
instability, comparable to the case of standard hybrid in- 
flation. The end of inflation tp cn d hence becomes a free, 
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FIG. 2: Evolution of the scalar field tp (solid line) and of the 
Lorentz factor 7 (dotted and dashed lines) as a function of 
the number of e-folds N, for a ~ 0.04 (m = 0.01m pl , A = 
10 ). The initial value of the inflaton field is such that 71m = 
3. The solid line represents the trajectory of the inflaton 
field computed from the exact equation of motion. The DBI 
slow-roll trajectory of Eq. (|111|> cannot distinguished from the 
exact one on this plot. Concerning 7, except during a short 
transient regime at the beginning of inflation, the two curves 
also match. 



additional parameter of the model, which is notably cru- 
cial for normalizing predictions to CMB observations [5J| . 
The accuracy of the DBI slow-roll trajectory can be 
assessed by an exact numerical integration of the equa- 
tions of motion. The results are plotted in Fig. [2] and 
confirm that, except during a few e-folds after the be- 
ginning of inflation, the DBI slow-roll approximation is 
accurate. Similarly, the other Hubble flow functions can 
be approximated by 
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where Eq. I|19]) has been used. As a result 
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Concerning the sound flow functions, one gets 

-3/2 



27TQ! \ m 4 . ' I 



S 1 ~ — M- + - 



3_ 

27T V rn 



x 4 / 4 

4 



111 ' 



-3/2 



Similarly, for small values of tp 
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FIG. 3: Exact numerical evolution of the Lorentz factor 7 
(solid line) and ei as a function of the renormalized field value 
a ' (p/m pl and for various values of the parameter a. Notice 
that inflation proceeds from larger towards smaller field val- 
ues. As can be seen from Eq. (|114p . the a-dependence of 
the Lorentz factor 7 is in the renormalized field. One can 
check that the asymptotic values of ei are compatible with 
the slow-roll predictions given in Eq. (|115p (see also Fig. U}. 
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FIG. 4: Exact numerical evolution of the slow-roll param- 
eters as a function of the number of e-folds N. The other 
parameters are the same as in Fig. [3] 



inflationary flow formalism (besides the fact that it is 
hard to interpret since there is no natural measure in pa- 
rameter space) where the Hubble scale H(<p) and 7(95) 
are expanded in terms of the inflaton field at a specific 
order |28j . 

Owing to the above analytical expressions, we can now 
easily consider the observational predictions of this sce- 
nario. The multipole moments at large angular scales are 
given by 
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which, for i = 2, gives the quadrupole moment 
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The various quantities appearing in the above expressions 
should be evaluated at "sound horizon crossing" . Using 
the trajectory (|113p 
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where N* is the number of e-folds between the end of 
inflation and the sound horizon crossing. It is usually 
considered that 40 < N* < 60. However, for the model 
under consideration, we have 7 ex l/<p 2 and H 2 ex ip 2 and, 
therefore, these two terms cancel out of Eq. (|124| . Since 
ei ~ y/a/(4n), we have in fact the remarkable property 
that the multipole moments do not depend on the num- 
ber of e-folds between sound horizon crossing and the end 
of inflation. This property was first noticed in Ref. [501 ] - 
Working out the previous expressions, one arrives at 
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Since a « 1, this means that the mass m should be 
smaller than in the standard chaotic scenario. The pre- 
vious expression can be re-expressed in terms of the pa- 
rameter A, 
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(128) 



Inserting Eqs. (|120[) and (|123|) into the expression of 
the scalar spectral index derived in Sec. Ill E H one ob- 
tains 



The evolution of the DBI slow-roll parameters is repre- 
sented in Fig. [3] and Fig. [4] and compared to an exact 
numerical integration of the equations of motion. Again, 
apart during the few e-folds of the transient regime, 
the DBI slow-roll trajectory and parameters are in good 
agreement with the exact integration. 

As already noticed before and as clearly shown in 
Figs. [5] and [3] the Lorentz factor 7 goes to infinity while 
<p — * 0. In this case, one may question the use of the 



l = e>(e 3 ,r5 3 ,e l <F), i+j 



(129) 



The spectral index vanishes at first order in agreement 
with the results of Refs. [20, [23] ■ Here, we find that, for 
this model, this also the case at second order. Notice also 
that from Eq. (|100|) . the second and third-order contri- 
bution in the running of the spectral index vanishes for 
if — » and a s reads 



0(e 4 ,<5 4 ,e l <P' 
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Let us notice however that, although extremely small, 
the actual value of a s obtained from the DBI slow-roll 
approximation is significantly different than the value ob- 
tained from the numerical integration, typically at almost 
100%. Such a loss of accuracy is related to the extremely 
flat power spectrum which has an almost vanishing run- 
ning; as for de Sitter, such a limit is in fact singular for 
the scalar modes. In any case, such a running is by far too 
small to be detectable in any of the present and planned 
cosmological experiments. 

Then, one can compare the above predictions with the 
WMAP5 0, 0, H H 11 Hi constraints on the CMB power 
spectra established in Ref. 45J. Notice that, in order to 
perform this comparison, it is mandatory to use a theo- 
retical power spectrum which takes into account the fact 
that the sound speed is a time-dependent quantity. It 
would be inconsistent to do this comparison using the 
constraints on the ordinary slow-roll parameters estab- 
lished, for instance, in Rcfs. [5J, (55|. From Ref. [4a |. 
we see that, in absence of a significant contribution of 
gravitational waves, a red tilt is favored. Therefore, the 
present model which is scale-invariant to a very high ac- 
curacy is clearly disfavored by the current data. 

Let us remind that the model is also strongly disfa- 
vored by the WMAP non-gaussianity bounds, as shown 
in Ref. [20| . Indeed, at lea ding order for DBI models, the 

Hi 



parameter / NL is given by [26 
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which, in the present context, gives 
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because M > 1 and a < 1. The above inequality is in 
contradiction with the observational range —151 < / NL < 
253 (see Ref. @). 



3. Chaotic Klebanov Strassler models with constant term 

We now consider the case where the Coulomb term is 
negligible and the potential can be written as 



V{ip) = V Q + -m\ 2 , 



(134) 



while the warp function remains unchanged and given 
by its expression in Eq. ([55]) . These models have been 
discussed in Refs. 0, S3, 11 H HH . 

Using Eqs. dUD and Q1J) leads to <p = -sm 2 p/(3jH). 
This means that, in the UV (s = 1) case, the vev of 
the inflaton field is decreasing as inflation proceeds while 
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FIG. 5: Evolution of the Lorentz factor 7 and of the Hubble 
and sound flow functions in the IR case of the CKS models 
with constant term. The parameters are a = 500 and f3 — 
1. Inflation proceeds from left to right and the "divergence" 
occurs at ip/m pl — > \/2/3. 



it is increasing in the IR case (e = — 1). In order to 
understand in which regimes DBI inflation can occur, it 
is convenient to express the Lorentz factor in terms of 
the vacuum expectation value of the inflaton field as 
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(135) 
Of course, when e = +1 and = 0, one checks that 
Eq. (|114[) is recovered. The DBI regime, 7 — > +00, is ob- 
tained when ip — ► 0. In this case, 7 can be approximately 
written as 
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Let us notice that this last equation is valid for e = ±1. 
In addition, in the IR case, the Lorentz factor also blows 
up when ip — * m P1 -\/2/3" as 
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The evolution of 7(cp) is displayed in Figs. [5] and [5] No- 
tice that we must always have <p < m pl \f2J3 in order to 
guarantee the positivity of the potential. Moreover, since 
in the limit ip — ► rripjV^S, V(ip) —> 0, the corrections to 
the potential may become non-negligible. In particular, 
the Coulomb part of the potential that we have neglected 
in this subsection will become important again. Notice 
that this situation is rather similar to the case of small 
field models in standard inflation where the potential is 
given by V oc 1 — </> 2 and where, at the end of the slow- 
roll phase, it is necessary to add higher order terms of 
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the form (j? in order for the potential to have a mini- 
mum which ensures that the reheating can proceed. We 
conclude that the limit <p — v m Pl -\/2j3, although it does 
imply 7 — > +00, appears to be more difficult to realize 
from the physical point of view. 

For all e, the first slow-roll parameter can be written 
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which, in the limit ip — > 0, becomes 




(139) 



and goes to zero as the inflaton vev vanishes. This means 
that, in the UV case, where the field decreases from some 
initial value towards zero, inflation cannot stop by viola- 
tion of the slow-roll condition. 

On the contrary, in the IR case, the field starts from a 
value close to zero and increases. As already mentioned, 
the vev of the field is bounded by \f2j3. In this limit, the 
first slow-roll parameter blows up as 
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and contrary to the UV case, inflation will stop by vio- 
lation of the slow-roll conditions. The evolution of ei is 
displayed in Fig. \E\ (IR case) and in Fig. [5] (UV case). 

It is also interesting to establish the expression of the 
other slow-roll parameters. Lengthy but straightforward 
calculations lead to the following expressions 
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These functions are represented in Fig. [5] for the IR case 
(e = —1) and in Fig. [S]in the UV case (e = +1). 



With a constant term included, Eq. ([20]) describing the 
field trajectory becomes significantly more complicated. 



10' 

£ 10 1 

CO 

g. ict 1 



o 10" 

CO 

^ 10" 
10" 



UV 



0.0 



0.1 



0.2 0.3 

¥>/m pl 



0.4 



0.5 



FIG. 6: Same as Fig.[5]but in the UV case. The inflaton field 
decreases while inflation is under way. 



Expressed in terms of the parameters a and (3, the DBI 
trajectory reads 



N(<p) = - 2ne/3 / du(u- e) 3/2 
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where we have defined the new dimcnsionless variable u 
by 
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(146) 



The expression (|145[) is an elliptic integral and can be 
expressed in terms of the canonical elliptic integrals of the 
first and second kind, respectively denoted as in Ref. [57| 
by F(9 \p) and E(8 \p). We will further use the definition 
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To ensure that all expressions are well-defined in the fol- 
lowing, we now consider separately the case of IR and 
UV models. 

For the IR case, e = — 1 and it follows from Eq. (|146p 
that the variable u maps the range p/m p] £ [0, \f2ft\ to 
u e [-1, 1]. From Eq. (| 147ft . we see that < n < 1 and 
both u and n appearing in Eq. (|145|) are in the canonical 
domain of definition of the elliptic integrals E(8\p) and 
F(6 \p) (see Ref. J57J]). One gets 
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FIG. 7: Typical evolution of the scalar field <p according to 
Eq. (|148[1 in the IR case. The parameters are a = 4 and j3 = 
0.05 such that the condition (3/a <^C 1 is satisfied. The initial 
condition is <p ini /rn pl ~ 0.022 and one can check that the field 
value always remains sub-Planckian. With the parameters 
considered here, one has n = 0.1, i.e. < n < 1 as required 
for the IR case. 



A typical IR trajectory is represented in Fig. \7\ This 
expression, although exact, is quite involved. Therefore, 
it is interesting to consider the dominant behavior for 
u — > —1, or <p — > 0, i.e. for very early times. In this situ- 
ation, it turns out that the terms containing the square 
roots are the dominant ones because they contain a pole. 
Then, neglecting the other terms, one can express the 
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variable u and, hence, the inflaton field vev p in terms of 
the number of c- folds. This leads to 
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One checks that when the number of e- folds increases, the 
vev of the inflaton field increases. This clearly reproduces 
the behavior seen in Fig. [7J Notice that, in order to 
establish the previous formula, we have also assumed that 
Pini/m-p-i <C 1 since p/m m <C 1 and (p- m { < (p. Moreover, 
the approximation is valid only when the second term in 
the squared bracket of Eq. (|149|) is small in comparison 
to one such that the vev remains positive. When N is too 
large, the brane is far from the bottom of the throat and 
this approximation breaks down. Finally, in the previous 
section, we have established the constraint (3/a <C 1. We 
see that, in Eq. (|149p . the inverse of this factor appears 
and in order to be consistent with the condition ip/m m <ti 
1, one must also have ipi n i/m pi <C y/0/a. 

Let us now turn to the UV case. In this situation, the 
variable u covers the range u6 [1, +oo[ where the lower 
bound corresponds to <p — ► 0. This is out of the canonical 
domain of definition for the elliptic integrals and to carry 
out the integration of Eq. (|145[) . one can introduce the 
variable 1/u. Moreover, the parameter n has now two 
disjoint ranges: for < a/3 2 < 1, n s] — oo,0[, whereas 
for a(3 2 > 1 we get n g]l,+oo[. Some care is therefore 
required on the formulations of the elliptic integrals. 

We first focus on the case where a(3 2 > 1. After ap- 
propriate redefinitions [57|, the final trajectory reads 
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At the end of inflation, i.e. when ip — > (or u — > 1), the 
next term to last in Eq. (|150p dominates, and we can 



invert the resulting expression to obtain 
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Let us notice that, this time, the factor \/(3/a <C 1 di- 
rectly appears in the above expression and, therefore, 
guarantees the consistency of the formula expressing the 
quantity tp/m pi <C 1. 

Let us now consider the case < a/3 2 < 1. The ellip- 



tic integrals F(6 \p) , E{9 \p) that arise in the calculation 
of Eq. (|145|) have now to be redefined for complex ar- 
guments 9 and/or negative parameters p, respectively. 
After some calculations [57|, the trajectory reads 
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In the limit u — > 1 or <p — > 0, the next to last term in 
Eq. (|152p dominates and we can approximately invert the 
above expression. The resulting trajectory ends up being 
given by the same expression as in Eq. (]151|1 . Finally, two 
typical UV trajectories are represented in Fig. [5] 

Using these results, let us now turn to the observational 
predictions of these models. The Cosmic Background 
Explorer (COBE) normalization, implemented with the 
help of Eq. (|125p . leads to the following relation, valid 
for the UV and IR cases 
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This formula is only consistent when (p*/m pl is a small 
quantity. In the limit ip — > 0, one has from Eq. Q139P 
ei oc ip 3 and 
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Therefore, e 2 and 5\ give the dominant contribution to 
the spectral index and one has 

n s - 1 ~ 4<?i. (155) 

Using the COBE normalization given above leads to 
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The spectral index only depends on the dimensionless 
't Hooft coupling constant A which is quite remarkable 
given the complexity of the equations and the number of 
free parameters of the model. In the UV case, the tilt is 
positive while it is negative in the IR case. In order to 
compatible with the CMB data, one sees that one must 
have A > 3 x 10 13 . 

The calculation of the running proceeds along the same 
lines. In the limit if — > 0, one has 

S 2 ~ -6i, 6 3 ~ \(l-4e)8i. (157) 

Then, using Eq. (fTTj) . one obtains the following expres- 
sion 
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The IR models (e — —1) have therefore a red spectral 
index and a positive running while, on the contrary, the 
UV models (e — 1) have a blue spectral index and a neg- 
ative running. In Fig. [51 we have plotted the scalar power 
spectrum for the UV model stemming from an exact nu- 
merical integration of Eq. (|24| . Both the spectral index 
and the running match at a few percent with Eqs. (|156p 
and 1(158]) . 

One can also estimate the tensor to scalar contribu- 
tion. For this purpose, one can compute the parameter 
r = 16ei/7, see Eq. (JT3"]) . Using again the COBE nor- 
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FIG. 8: Typical evolution of the scalar field ip according to 
Eq. (I150|l and Eq. (|152|l in the UV cases. The solid line refers 
the UV case with n > 1 since the parameters chosen are a — 
404 and /3 = 0.05 which implies n = 101. On the contrary, 
the dotted line represents a trajectory in the UV case with 
n < 0. The corresponding parameters are a — 4, P = 0.05 
which gives n = —0.01. In both cases, the initial condition is 
(/3 ini /m pl ~ 0.67. 
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FIG. 9: Exact (solid) and first order (dotted) scalar power 
spectrum for the CKS plus constant UV models. The parame- 
ters are a = 38, f3 = 3.7, A = 1.6xl0 13 and m 2 = SxlO" 4 ™.^. 
From the exact spectrum, one gets at the pivot n s — 1 ~ 0.11 
and a a ~ —0.0023 in agreement with Eqs. (|156|l and (|158|l . 



IV. CONCLUSION 

In this section, we briefly recap and discuss the results 
obtained before. The main goal of our article was to 
present a generic formalism, as similar as possible to the 
standard approach, to investigate the features and the 
observational predictions of k-inflation. In particular, we 
have offered a new formula (but valid for DBI inflation 
only) for computing the classical trajectory for any warp 
function T(ip) and potential V((p). Then, using the uni- 
form approximation, we have computed the scalar and 
tensor power spectra, see Eqs. ([5U)) and ([57)1 . These for- 
mulas are based on a double expansion. There is first an 
expansion around a pivot scale and, then, each coefficient 
of this expansion is in turn evaluated by means of another 
expansion, this time in terms of the Hubble and sound 
flow parameters. Finally, we have illustrated our results 
on various example models encountered in the literature. 
In particular, we have derived new results concerning the 
CKS models with constant term and obtained simple ex- 
pressions for the cosmological observables. 

Having generalized the standard formalism to k- 
inflation, the next step is of course to compare the pre- 
dictions to the observations. This is done in Ref. [45| in 
the case of the CMB WMAP5 data. 
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APPENDIX A: POTENTIAL AND WARP 
SLOW-ROLL PARAMETERS 



In a broader sense, one can also try to establish an- 
other hierarchy of expansion parameters to the DBI case, 
namely in terms of the potential V(tp) and the geometry 
parameter T((p). If we define the two parameters e v and 
e T according to 



malization, straightforward manipulations lead to 
_ 128^T c 2 mb ! 



15Q 2 



(159) 



which is very small since f3/a <§; 1 and A 1 ex (n s — l) 4 . 
Finally, the level of non-gaussianity can be estimated to 



JNL — ao a ' 



/45 Q m pl /m 



-9.2 x 10 



which is negative. 
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(160) 
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where, in the present context, a prime means a derivative 
with respect to <p. We can establish the link with the 
(ej, Si) parameters. Straightforward manipulations lead 
to the following expressions 



7 SI 
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l + 7(l-ei) 



(A2) 



(A3) 
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As can be expected, there is no one-to-one correspon- 
dence between the two sets of parameters (i.e. , for in- 
stance, the expression of e T not only involves ei but also 



APPENDIX B: RELATIONS WITH OTHER 
SLOW-ROLL HIERARCHIES 



the constant C in that reference being related to the con- 
stant D used in the present paper by D = (C — 3)/4. 
The small difference in the actual numerical values of 
(C - 3)/4 = In 2 + 7Euior - 7/4 ~ -0.73 and our 
D = 1/3 — In 3 ~ —0.76 is the just the well-known dif- 
ference between the WKB approximation versus the first 
order integration of the Bessel equation [43[ ■ 



In this Appendix, we briefly state the relations between 
our hierarchy of (ej,<5j) parameters and the various pa- 
rameter sets that have been used elsewhere in the lit- 
erature. The parameters (ec,?7_D, «D)£dj Pd) defined in 
Ref. [35| can be expressed in terms of our parameters as 
follows: 



en = 


= ei, 


(Bl) 


Vd -- 


= ei - -(e 2 + <5i), 


(B2) 


n D = 


= Si, 


(B3) 
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£2£3 . Sie 2 3 

= ei £2 
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+ei f ei - - Si J + — (<$i + <5 2 ), 


(B4) 


PD = 


= -J-(2& + 3<Ji-e a ). 


(B5) 



In Ref. [34| , the first three of these parameters were used, 
changing the notation to to — > e and so on. Then, it is 
interesting to compare the expression (2.30) of Ref. [35[ 
to our Eq. (|58"j). Written in terms of the Hubble and 
sound horizon parameters, this equation reads 

n s - 1 = -2e l0 - e 2o + S l0 - 2ef + (2D - 3)ei e 2o 

+3eio^io + £2o<5io + -De2©e3o - <5i 
+2DS 1o S 2<> . (B6) 

The claim made before can now be explicitly checked. 
The above formula coincides at first order with Eq. ((68]) 
but differs at second order. It is interesting to notice that 
the differences show up only in those terms containing the 
constant D. 

A slightly different set of parameters (e, s, rj, p, 2 A) was 
recently introduced in Ref. [27|, and expressed in terms 
of our (ei,Si) these parameters read 



e = 


ei, 


V = 


ei - 


s = 


Si 


p = 


Si 
2ei 
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eU 



;{e2 + Si), 



(2S 2 + 3Si-e 2 ), 

7:(e2e3 — 3eie 2 + e 2 5\ - 3ei^i 



+Sf + S t S 2 ). 
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(BIO) 

(Bll) 



Using this correspondance, one can explicitly check that 
Eq. (84) of Ref. [27| matches exactly with our Eq. 



APPENDIX C: COSMOLOGICAL 
PERTURBATIONS FROM DBI INFLATION 

In this Appendix, we briefly establish the conservation 
law for cosmological perturbations in k-inflation. This 
is important because as in the standard case 58j, this 
allows us to propagate the primordial power spectra to 
the recombination epoch. We also briefly recall how the 
equation of motion for the Mukhanov-Sasaki variable can 
be obtained. 

In the longitudinal gauge, the perturbed FLRW metric 
reads 



ds 2 = a 2 (r;){-(l + 2$)d7 ? 2 

+ [(l-2^)6 ij +h ij ]dx i dx j }, 



(CI) 



where $ and \t are the Bardeen potentials which are 
coupled through the perturbed Einstein equation to the 
field perturbations Sip. The transverse and traceless ten- 
sor hij, i.e. satisfying hi — d°hij — represents the spin 
two fluctuations which are not considered in the following 
since they are not affected by the DBI dynamics. 

At the first order, the perturbed energy momentum 
tensor is obtained from Eq. ([1]) and its components read 
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(C4) 



The gravitational sector remaining standard General Rel- 
ativity whose perturbed Einstein tensor can be found in 
Ref. [59|. In particular, one still has the relation ^ = $ 
and there is only one degree of freedom because $ and 
Sip are related by the perturbed Einstein equations. One 
can therefore reduce the study of the scalar sector to the 
study of a single variable, as the comoving curvature per- 
turbation 



c = $ 



H 



n 2 -n 
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(C5) 



From the background and perturbed equations, one can 
show that the comoving curvature perturbation can be 
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simplified to 



( = <P + H 



Sep 



(C6) 



which also matches with its usual expression in standard 
inflation. Straightforward manipulations allows us to de- 
rive the expression of the derivative of C and one gets 



c 



H 



H 2 - H' 7 2 



$. 



(C7) 



As a result, £ is a conserved quantity on scales larger 
than the sonic horizon and allows us to propagate the 
spectrum from horizon exit till the beginning of the radi- 
ation dominated era. As usual, this result applies if the 



decaying mode is neglected and in absence of entropy 
(isocurvature) perturbations [58| . 

The Mukhanov-Sasaki variable is now related to the 
comoving curvature perturbations by «k = zQ where 



7' 



3/2 a V 



n ' 



(C8) 



has a 7 dependence. Using the previous equations, Vk is 
found to obeys the mode equation 



k 2 z" , 

-2 Vk = 0. 

7 Z z 



(C9) 



This is the equation already derived in Refs. [27], [34J and 
used in the text with the correspondence c s = I/7. 
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